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Abstract 

We evaluate charge diffusion constant of dense and hot hadronic 
matter based on the molecular dynamical method by using a hadronic 
collision generator which describes nuclear collisions at energies 10 1 ~ 2 
GeV/A and satisfies detailed balance at low temperatures (T < 200 
MeV). For the hot and dense hadronic matter of the temperature 
range, T = 80 ~ 200 MeV and baryon number density, ns = 0.16 
fm~ 3 ~ 0.32 fm~ 3 , charge diffusion constant D gradually increases 
from 0.5 fmc to 2 fmc with temperature and is almost independent 
of baryon number density. Based on the obtained diffusion constant 
we make simple discussions on the diffusion of charge fluctuation in 
ultrarelativistic nuclear collisions. 
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1 Introduction 



Searching for the quark gluon plasma (QGP) is one of the hottest topics in 
the recent high energy nuclear /particle physics |J. Many kinds of phenom- 
ena have been proposed as candidates for the experimental signal to find 
the QGP. Charge fluctuation belongs to one of the most promising signals. 
Occurrence of Disoriented Chiral Condensation (DCC) at QCD phase transi- 
tion would cause the large fluctuations of the ratio of the numbers of charged 
pions and neutral pions [TJJ. Difference of the fluctuation intensity between 
QGP phase and hadronic phase can work as a signal of the existence of QGP 
H f|, [5|]. However, even if QCD phase transition takes place and the char- 
acteristic fluctuations of QGP or QCD phase transition are produced, such 
a fluctuation can be wiped out during hadronic era before the freeze-out. 
Whether the fluctuation which is caused at the phase transition or produced 
in the high temperature phase can survive in the hadronic era or not is the 
diffusion problem or the transport problem in the hot and dense hadronic 
matter. Macroscopic phenomenological equations, e.g., Navier-Stokes equa- 
tion and diffusion equation, enable us to describe such phenomena in a simple 
manner. However, those kinds of phenomenological equations contain mate- 
rial constants, so called transport coefficients, and the dynamical calculation 
of the transport coefficients is a very important and difficult problem of sta- 
tistical mechanics. 

Interactions between hadrons are strong interaction which are believed to 
be described by QCD. However, in the hadronic energy region, perturbation 
does not work for QCD and no systematic way to treat is established. In a 
previous paper , we have reported a calculation of baryon number diffusion 
constant based on the relativistic collision event generator URASiMA (Ultra- 
Relativistic AA collision Simulator based on Multiple scattering Algorithm) 
which can reproduce hadronic spectra of nuclear collisions in the BNL-AGS 
and CERN-SPS. Usually collision event generators are so designed as to be 
suitable for the description of multiple production, and the detailed balance 
between the interactions is used to be paid only little attention [7). As we 
discussed in the ref. ||, we improved the URASiMA to recover detailed bal- 
ance in the hadronic time scale and we have succeeded to establish stationary 
states of interacting hadrons with fixed temperatures and fixed baryon num- 
ber densities ||. In this paper we evaluate charge diffusion constant in the 
hot and dense hadronic matter and investigate the diffusion of the charge 
fluctuation in the relativistic nuclear collisions. 
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2 Statistical ensembles 



In order to prepare the statistical ensembles for the hot and dense hadronic 
state, we put hadrons in the box and updated with numerical code URASiMA 
with periodic boundary conditions. Recipes for the simulation are quite 
similar to the ordinary molecular dynamics. Molecular dynamics is one of 
the well established numerical methods in the statistical physics, however, it 
has been developed mainly for nonrelativistic system where particle number 
is conserved. The hot and dense hadronic system in which we have interest 
is a fully relativistic system, and the particle production and decay occur 
naturally. Therefore, the existence of stationary state itself is not apparent 
in this system. 

URASiMA is originally designed as an event generator for the relativis- 
tic nuclear collisions based on the hadronic multichain model (MCM). In the 
high energy collision experiment, multiple production takes place and the sys- 
tem is thought to expand quickly. Therefore, usually, production processes 
play essential roles in the relativistic event generator but the reabsorption 
processes (reversal process of multiple production process) do not have been 
thought to be important, and sometimes have been neglected in the simu- 
lation code. However detailed balance of interactions is very important for 
statistical physics and naive application of collision event generator to the 
molecular dynamics in the box leads to the one-way conversion of the energy 
into particle production. As a result, Hagedorn type behavior appears, i.e., 
strange saturation of the temperature occurs ||, [l(| . 



In order to recover the detailed balance, we have improved URASiMA 
to contain many resonances and changed the code so as to describe some of 
production processes occur through production and decay of the resonance. 
The reversal processes of those processes have been naturally taken into ac- 
count. As a result, after initial thermalization time period of about 150 fm/c, 
detailed balances among interactions seem to be almost kept during the sim- 
ulation. It is noted that such improvement does not spoil the descriptive 
power for ultrarelativistic nuclear collisions ||. As we have already reported 
in ref. @, slope parameter T of the energy distribution (Table [I]), 

dN dN 

d^ = 4^dE =CeM - E/n 

became almost common value for all particles (fig. 1) and the population of 
particles became stationary (fig. 2). Therefore, we looked upon the system as 
the equilibrium state with temperature T and fixed baryon number density. 
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Table 1: Fitted parameter T of each particles at t = 150 fm/c. Baryon 
number density is normal nuclear density, ub = 0.156 fm -3 . 



e tot [GeV/fm 3 ] 


N 938 [MeV] 


Ai 232 [MeV] 


7i [MeV] 


P770 [MeV] 


0.313 


131 ± 5 


122 ± 3 


132± 2 


141 ± 5 


0.625 


170 ± 6 


159 ± 5 


163 ± 1 


168 ± 2 


0.938 


181 ± 6 


177 ± 6 


187 ± 1 


190 ± 2 



fig. 1 



fig- 2 



Running URASiMA many times with the same energy and the same 
number of baryons, Nb, in the box with volume V, and taking equilibrium 
configurations, we have prepared statistical ensembles of the state with tem- 
perature T and baryon number density Nb/V. Throughout this paper, we 
assumed iso-symmetry; the same number of protons and neutrons are put at 
initial time. 



3 Diffusion constant in the linear response 

According to Kubo's Linear Response Theory, diffusion constant D is ob- 
tained current (velocity) correlation Because of the relativistic property 
of the hot and dense hadronic state, we should use f3 = p/E instead of usual 
v, 

1 roo 

D = - < (3(t) ■ (3(t + t') >dt'c 2 , (1) 
3 Jo 

with c being the velocity of light. In the calculation of the charge diffusion, 
average is taken over all charged particles. When we evaluated baryon num- 
ber diffusion constant Db in our previous paper, we took average over only 
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baryons. If correlation of velocities damps exponentially, 

t' 

</3(t)-/3(t + 0>ocexp( — ) (2) 

T 

with r being relaxation time, we can rewrite eq. (|l|) with the simple form as, 
D = i < (3(t) ■ f3(t) > c\ (3) 

with use of the relaxation time. 



fig. 3 



Figure 3 shows the correlations of the velocity of charged particles and ex- 
ponential damping seems very good approximation. Diffusion constant ob- 
tained through eq. (EH is shown in fig. 4. 



fig. 4 



Diffusion constant increases with temperature and is almost independent of 
baryon number density. This result shows clear contrast against baryon num- 
ber diffusion in our previous paper. Baryon number diffusion constant shows 
baryon number dependence and changes with temperature only weakly. Both 
calculations are almost the same but main contribution to each transport is 
different and we can understand fig. 4 as follows: In charge diffusion, main 
contribution comes from charged pions of which mass is just comparable to 
the temperature of the systems. On the other hand, the mass of the baryon 
is much larger than the temperature and baryonic system is almost nonrel- 
ativistic. The number of the baryon (baryonic charge carrier) is determined 
by the baryon number of the system, and even in very low temperature re- 
gion, baryons must exist because of the fixed baryon number. In the lower 
temperature than pion mass, pion degrees of freedoms are frozen and elec- 
tromagnetic charge and baryonic charge are both carried by baryons. At 
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temperature about pion mass, pion degrees of freedom start to be melted and 
start to contribute to the charge transfer. In the higher temperatures, pions 
dominate the charge transfer. As a result, in the temperature lower than 
pion mass, both diffusion constants are almost the same, but in the temper- 
ature higher than pion mass, diffusion constant of charge current gradually 
increases with temperature. 



Charge fluctuation in ultrarelativistic nu- 
clear collisions 



Using the above obtained values, let us sketch the diffusion of the charge 
fluctuation in the relativistic nuclear collisions. In the simplest static picture, 
the solution of diffusion equation, 

^f(x,t) = DV 2 f(x,t), (5) 




1 \ 3 (X-Xg) 2 

e 2R 'o (6) 



2itRl 



3 jX-X a f 

e 2fl2+4D(i-i ) _ (7) 



n(2R 2 + AD(t-t )) 
Suppose that charge fluctuation with size Rq produced at t — to, subjecting 



to eq. (^|), the fluctuation diffuses and expands to R(t) = \JRq + 2D(t — t ). 
Therefore, we may regard that at such time t that 2D(t — to) — Rqi charge 
fluctuation almost disappears. If the charge fluctuation about Ro = 3 fm 
at T = Tq = 160 MeV is produced, it survives only 2 fm/c during hadronic 
matter era. But if charge fluctuation is produced at lower temperature about 
120 MeV (Super Cooled DCC), diffusion constant is about 0.7 fmc and the 
fluctuation with initial size of about 3 fm will survive for 6 fm/c. The length 
of hadronic era and the chronological change of the temperature depend on 
the solution of the hydrodynamical model, but the latter case is promising 
to be observed experimentally. 
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In order to take account of expansion of the hadronic fluid, we must use 
relativistic hydrodynamical model. In the relativistic hydrodynamical model, 
charged current is given as, 

J* = f{x,t)U» + J£, (8) 



with being local four velocity |13|. f(x, t) is charge density on the local 
rest frame, f(x,t) = U^J^. The diffusion current, J% is given by relativistic 
extension of the Fick's lawQ 

J2f = -AW(x,i) I (9) 
c 

where A^ v = g^ v — U fJ, U u is space-like projection operator orthogonal to £/ M . 
Putting above J M into the continuity equation, d^J^ 1 = 0, we can obtain the 
relativistic diffusion equation, 

U»dJ(x, t) + f(x, t)dJJ» + ^d^ u d v f(x, t) = 0. (10) 



In the case of 1+1 dimensional scaling expansion, particular solution of 
the local four velocity is given as C/ M = — with light-like coordinate, ct = 
ct cosh r), z = ct sinhr/ JTJ]], and diffusion equation (|I0|) is rewritten as, 

§-f(v,r) + lf(v,r)^^^- 2 f( V ,r). (11) 



The solution of eq. ( |TT| ) with initial condition, f(rj, r ) = y^^z- exp(— ^ )/ , 
is given by 

'('-^tSf^'' (12) 

2 



with -R(r) being 0, 



R(r) 2 = Rl + 2^(r-r ). (13) 



The first factor in eq. (|12~D , -, is the result of systematic expansion of the 
scaling solution; the region with fixed At] expands as tAt]. This rapid expan- 
sion also dominates the evolution of the fluctuation size, R{t) 2 ; the region 
with fixed Arj expands with r but the diffusion effect subject to eqs. (||) or ( 

-"in the relativistic notation, x° = ct, and = -M^. 

2 Note that this is the difussion in the rj direction and R is the size in the rj coordinate. 
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TUP is the expansion with only y2Dr. Therefore, in rj coordinate, the effect 
of the diffusion becomes smaller and smaller with time, According to 

our result (fig. 4), diffusion constant in the hadronic matter is at most 2 
fmc and the diffusion term in eq. ( PH) becomes small enough already at r = 
several fm. As a result, the fluctuation in rj will be deformed only little by 
the diffusion effect in the hadronic era. 



5 Concluding remarks 



We evaluated charge diffusion constant of hot and dense hadronic mat- 
ter based on the relativistic collision event generator URASiMA. Obtained 
charge diffusion constants increase from 0.5 fmc to 2 fmc in the temperature 
range between 80 MeV and 200 MeV. Based on the obtained diffusion con- 
stant, we made rough sketches of the diffusion of charge fluctuation in the 
ultrarelativistic nuclear collisions. For more improved discussion, we need to 
solve the hydrodynamical equation coupled with the charge current conser- 
vation equation Jf5| . 
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n B = 0.156 [fm ],t = 150 [fm/c] 
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e = 0.31 3 GeV/fm 3 e = 0.625 GeV/fm 3 



e = 0.938 GeV/fm 
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Figure 1: Energy distributions of particles at t — 150 fm/c. Lines stand for 
the results of fitting with thermal distribution. Normalizations of the data 
are arbitrary. 
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Figure 2: Particle density of time. The population of the particles in later 
than t = 150 fm/c seem to be stationary. In each case, baryon number 
density is n# = 0.156 fm~ 3 , and energy density for (a) is 0.313 GeV/fm 3 , (b) 
is 0.625 GeV/fm 3 and (c) is 0.938 GeV/fm 3 , respectively. 
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Figure 3: Velocity correlation of the charged particles as a function of time. 
Lines correspond to the fitted results by exponential function eq. (2). Nor- 
malizations of the data are arbitrary. 
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Figure 4: Diffusion constant D and baryon diffusion constant D B - 
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